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ABSTRACT 


Higher phase order transmission (HPOT) is an 
important area of research in power systems. This is 
mainly due to the increasing demands for power which 
have forced engineering to explore newer avenues of 
power generation and transmission. Design of such multiphas 
systems requires that adequate protection schemes 
be also provided for. The thesis presents one such 
scheme for protection of multiphase systems. The 
scheme is known as ultra high speed protection 
relaying* It is based on the concept of travelling 
waves which are generated when a fault occurs on a 
system. 

Rapid fault clearing has many advantages. 

Important among these is the improvement in transient 
stability of the system because the kinetic energy 
stored in the system varies as the square of the 
clearing time. In the past ultra high speed relaying 
was not feasible because of absence of fast acting 
circuit breakers. With advances in breaker 
technology, IKS relaying has become feasible. 

The thesis discusses a frequency domain model 
for six phase transmission line* This model is suited 
to study of travelling wave generation and propagation* 
Fault analysis of multiphase systems has been studied 
and studies conducted on a sample power system. An 



amplitude comparison relaying scheme based on 
travelling waves has been discussed along with its 


digital implementation* 



CHAPTER 1 


INTRODUCTION 

High phase transmission (HPOT) system is 
receiving increasing attention today as a response 
to the ever growing demands for power* Increase of 
power transmission by increasing system voltage seems 
to have reached a saturation point. High voltage 
has many problems associated with it. These include, 
corona discharge, radio interference and also 
biological hazards. But the most important reason 
which prompted research into multiphase systems is 
the difficulty in obtaining right of way. Due to 
increasing pressure on land and its Increasing costs, 
public utilities are finding it increasingly difficult 
to obtain new rights of way for transmission lines. 

It is hoped that development of multiphase systems 
would alleviate this problem to a large extent. 

Adequate protection schemes are essential for 
proper operation of any power system. Any protective 
scheme must satisfy the dual requirements of speed 
and reliability. In this context, ultra high speed 
relays have gained increased importance. However, in 
the past use of such relays was hindered because of 
the absence of ultra high speed circuit breakers. 

But with recent developments [1,2], realisation of 
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UHS protective relaying schemes has become imperative. 

Design of a relaying scheme requires that the 
system under study be accurately modelled. Ultra 
high speed relays discussed in the thesis use the 
concept of travelling waves generated during a fault. 

In other words, the relays are based on the transient 
waveforms generated after a fault. As a result, 
such relays have a distinct advantage over the 
conventional distance relays. Distance relays are 
based on calculation of impedance using the 
fundamental sinusoidal component of complex post 
fault waveform. This requires that, high frequency 
components and the dc offset be filtered out* This 
requires development of accurate analog/digital 
filters which results in time delays* On the other 
hand, travelling wave relays use the transient 
waveform as their input and hence, only the prefault 
fundamental needs to be filtered out* As a result, the 
relays are fast in operation. 

The thesis aims at developing a new ultra 
high speed relaying scheme for the protection of 
multiphase systems based on the concept of travelling 
waves. For this the following approach has been adopted, 
1* Mathematical modelling of a six phase transmission 
line for studying generation and propagation of 
travelling waves* 



2. Multi— node fault analysis of power systems 

using the bus admittance matrix taking into 
account parameters like fault resistance and 
fault initiation angle. 

3. Development of a relaying scheme for protection 
of a six phase transmission line using the 
travelling waves as relay input signals, 

4. Digital implementation of the proposed relaying 
s cherae « 

Digital simulation of faulted transmission 
lines, (3-phase) has been discussed by Johns and 
Agarwal [3], Kothari et al« [4] and Hump age et al* 
[5], Fault analysis techniques for three phase 
systems were discussed by Johns and Agarwal, 

A novel protection scheme based on travelling 
wave concept was first discussed by Vittins [6], 

Later work in this area was done by Chamia and 
Liberman [7], Matele [8] among others, Johns and 
Agarwal [9] also proposed a scheme which operates in 
conjunction with a carrier communication channel and 
In carrier blocking mode, 

A differential relay based on the concept of 
travelling waves along with its theory, sensitivity 
analysis etc* was proposed by Takagl [10] and Akimoto 
et al. [11], Digital protection schemes were studied 
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by Mann and Morrison [ 12,13]. Ramamoorthy and 
Verma [14] described a relaying principle in which 
a discriminant is used to determine the direction 
of the fault* 

UHS relaying scheme developed by ASEA, Sweden 
have been listed at the Bonneville Power Administra- 
tion (BPA) 500 k V system. 

This thesis in all contains four chapters, 
starting with the present chapter One on Introduction. 

The chapter Two develops a model for six phase 
transmission line suitable for transient analysis. 
Frequency domain Laplace transform has been used for 
this purpose. An equivalent Pi-model which Is 
suited to multi-node analysis has been derived. The 
chapter also presents techniques for simulating 
various types of faults on a six phase system using 
the superposition principle. The frequency domain 
bus admittance matrix has been used to study the 
fault generated voltages and currents, A sample 
3-bus six phase system has been used to test the 
techniques discussed in the chapter. The waveforms 
of fault generated voltages at relay location for 
various values of fault impedance has been presented. 

A new technique for Laplace transform inversion has also 
been presented. 

The Chapter 3 presents the proposed ultra high 
speed relaying scheme. The scheme uses amplitude 
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comparison of the two relay signals to decide whether 
the fault is internal or external* The relay input 
signals are derived from the modal voltages and 
currents* The chapter also discusses the digital 
implementation of the proposed relaying scheme* 
Finally the relaying scheme is simulated on a 
sample 3-bus power system. 

Finally Chapter 4 concludes with a review 
of the work presented in the thesis and scope for 
future work* 



CHAPTER 2 


DYNAMIC MODELLING AND TRANSIENT FAULT ANALYSIS OF 

MULTIPHASE SYSTEMS 

2.1 INTRODUCTION 

The chapter presents a model for a six phase 
transmission line which is suitable for transient 
analysis and travelling wave study. The model is 
developed in the Laplace S-domain. An algorithm 
which is used for inversion to time domain is also 
presented* The chapter also presents fault simulation 
of various types of faults on a six phase system and 
the algorithm for fault analysis of multinode systems. 

The model for a six phase transmission line 
developed is the S-domain it-model which is suitable 
for nodal analysis. In the past, frequency domain 
studies have been largely confined to the use of 
Fourier transforms. This was because Laplace transform 
inversion was cumbersome and the computations were 
unwieldy* The chapter discusses a new scheme for 
inversion of the Laplace transform. The algorithm 
discussed reduces Laplace inversion to a matrix multi- 
plication. The fault simulation invokes the 
superposition principle to calculate the fault 
generated components of voltages and currents • 
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To design a protection scheme for a given 
power system requires accurate determination of 
post fault voltages and currents. For this, accurate 
modelling of various power system elements is needed 
along with computational techniques which are 
numerically accurate as well as efficient. To this 
end, frequency domain analysis has been found to be 
more accurate [15], 

Modelling and fault simulation of three phase 
systems has been carried out by Humpage et al. [5], 
Johns and Agaxwal. However, these models suffer 
from the drawback that they are cumbersome and can 
be applied with ease only to simple systems. Modelling 
of multiphase power system elements has been done by 
Tiwari and Singh [16] and Chowdhary and Singh [17] • 
Most of these models are meant for steady state 
analysis such as load flow and planning studies etc. 

2,2 FREQUENCY DOMAIN PI-MODEL OF SIX PHASE 
TRANSMISSION LINE 

For simplicity in analysis, we assume that the 
line is perfectly transposed. This would imply, 

Z aa “ \b ~ Z ee " Z dd “ Z ee = z ff 

Z ab “ “ Z cd “ Z de “ *** “ Z fa’ 

Let the column vertex [V] denote the vector 
of the six phase to ground voltages, i.e. 
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V (s) 
a v 

V b (s) 

V c (s) 

V d (s) 

v e ( s > 

V f (s) 

Here, s is the Laplace domain parameter* Similarly, 
let [i] denote the column vector of six line currents 
in the six phases 

x a (s) 

I b (s) 

I c (s) 

I d (s) 

I e (s) 

I f («) 

Then, in the Laplace domain, the transmission line 
equations can be written as. 


d[V]/dx * 

- £z][i] 

U) 

d[j]/dx ** 

- [Y][V] 

(2) 


x is the transmission line's length parameter* 

Here [Z] is the impedance (per unit length) matrix 
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admittance (p.u« length) matrix, i.e. 

Z aa Z ab * 2 af 

Z ba ^bb 


[Z] « 



L Zfa 

and, similarly 






Y ab Y af 

Y bb 


[Y] 


'fa 


»»**#»*#*•■ 


l ff 


As stated above for a perfectly transposed line, we 
have , 


Z aa = Z bb = Z ec = ***** * z ff ** Z s 

Z ab = Z bc = Z cd = * Z fa * Z m* 

Using equations (1) and (2) above, we get 

d 2 [v]/dx 2 = [Z][Y][I] (3) 
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Equation (3) consists of six coupled ordinary 
differential equations. We should decouple these 
equations so that closed from solution can be obtained. 

It is easily seen, that, for a transposed line 
(assuming perfectly transposed line), the matrix, 

[P] * M [Y] 

is symmetric matrix whose off diagonal elements are 
equal and the diagonal elements are also equal* This 
suggests, that, we can use the linear transformation 
matrix (similar to symmetrical component transformation 
matrix for a 3-phase line) 


a=exp(— j60° ) 
s= 6f 1 

( 4 ) 

to decouple the above equations (3)* This trans- 
formation matrix [T] is linear, power invariant 
and complex* These decoupled equations are of the 
form. 
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T“ 

dx^ 


c V 

m 


(m) 


m = 0,1,2, 3,4,5 


( 5 ) 


is called the modal voltage and o m is 
the m tJfl mode propagation constant# Actually, 0 
represents the mth eigenvalue of the [p] matrix. 

Each of the equations (5) can be solved independently# 
Correspondingly, the modal current equations are. 


^2 (m) 

“"2 

dx 


= cr 1 

m 


(m) 


m = 0,1, 2, 3, 4,5 


solving equation (5) we get, 


(s) = K. cosh cr x + 1C, sinh x 

1 m z m 

m = 0,1, 2, 3,4, 5 (7) 

Then, using equation (l) we get, 

lW(s) ■ & (K 1 sinh v + *2 V) 

m « 0,1,2, 3, 4, 5 (8) 

Here, Z^ m ^ is the m^ h eigenvalue of the [Z] matrix# 
The boundary conditions at the receiving end are, 

1^ * 4 m) , and V (m) * V ^ at x * 1 • 

Similarly, at the sending end, 

l(m) = jW and y (m) = v (m) , t x _ 0> 

Using these boundary conditions, we get, for each 
modal component , 
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v<*>( S ) 


4 m) (s) 

33 

— ' 



(m) 


cosh cr X 
m 


m 


sinh cl 

m 


—T sinh cl cosh cl 


(m) 


m 


m 


L 


v| B) (s) 
4 b) (s)i 

( 9 ) 


Finally, the modal values are related to the 
corresponding phase values by the relation 


[V] 


[T] [V] 


m 


[I] - [T] [I] 


( 10 ) 


m 


[V] m and [l] m are the modal voltage and modal current 
column vectors respectively* 

Using the above transformation, we get, the 
following equations 

[V] H = [T][cosh cl][T]*[V] s - rntZ^rHsinh o£][T]*[I s ] 


■ [I] R = -[THoJtZjj]" 1 [sinh ni.][V s ] -[T][cosh <J«tT]*[I s ] 

( 11 ) 

where 

[V] R is the receiving end phase voltage vector | 

[V] s is the sending end phase voltage vector | 

[ij^] is the receiving end line current vector; 

[l s ] is the sending end line current vector; 

. [<3- ] diag[o Q c t a 2 a 3 c A o 5 J 

* [cosh oX] ** diag[coshe 0 X cosh cJL cosh a^l cosh o 3 X 

cosh c^l cosh 
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[sinh cr.t] = diag[sinh oJL sinh sinh GqQ. sinh cr ^ t 

sinh a^£ sinh a^.d ] 

and, 

[Z M ] = diag [Z^ 0 ^ Z^ Z^ Z^ Z^ Z^] 

The above equations can be written in compact form 


turn* p 



— 



•yps f Id p c p d $ 0 £ 

R 


[A*] 

[B*3 


v a,b,c,d,e,f 

V S 

▼a,b,c,d,e,f 

X R 


[C»] 

CD’] 


T a*b,c,d,e,f 

X S 


~ ( 12 ) 


The inverse equation can be written as. 


i/d,b ,c ,d, e ,f 
V S 


[A3 

[B] 


v a,b,c,d,e,f 

V R 

t a , b , c , d , e , f 
X S 

SZ 

[c] 

[D] 


«3| b , c , d , e , f 

X R A 

(1 


where 

[a] - [T]* [cosh oi] [T] - [D] 

(14) 

[B] - [T]* [Zy] [or 1 [Sinh o.t] [T] 

[C] - [X]* [o] [sinh oi] [ZyF 1 [T] 

The equation (13) represents the frequency S- 
domain equations which relate the receiving and 
sending end voltages and currents in Laplace domain* 
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Equation (13) forms the basis of the S-domain Pi-model 
of the six phase transmission line. 

Consider the Pi-model of the six phase 
transmission line as shown in the Fig. 1. Fran the 
model, we can write the following equations relating 
receiving and sending end quantities 


v a,b,c,d,e,f a(u + z a,b,c,d f e,f y a,b,c,d,e,f )v a,b,c,d,e,f 


'R 


. 7 a,b,c,d,e,f 7 a,b,c,d,e,f 

* 1 R 


I* a , b , c , d, e , f /^a,b,c ,d,e,f . *»a,b,c ,d,e,f 

V*2P + 


IP 


t y3jb jC jd| ©,f 7 a,b,c,d,e,f ^a,b,c,d,e,f \ 

’ 1 jjp ^*D *0 D * 


2p 


v a,b f c,d,e,f + (u + Y a.b f c,d,e.f 


7 a,b,c,d,e,f \ 7 a , b , c , d ^ e , f 
4 P ;I R 

(15) 


where, 

2 ra,b,c,d,e,f _ s ^ x p h ase transform impedance 

matrix of the series branch 
y a,b,c,d,e,f^ y a,b,c t d,e # f ^ six phase transform 

admittance matrix of 
the shunt branch 

U « Identity i.e, unit matrix of the order 6. 
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FIG. 1 FREQUENCY DOMAIN PI MODEL OF SIX PHASE 
TRANSMISSION LINE 





Comparing equations (13) and (15) and* we get 


c*d f e f £ 

y3|b|C 

^IP 

y3 |bf C 

X 2P 


= [B] 

- ([D] - U) [B]" 1 

= [B]' 1 ([A] - U) 


By expanding the matrices, it is easily seen, that, 

Y 1P = Y 2P* Thus we finall Y have the Parameters of 
the frequency domain Pi-model of a six phase trans- 
mission line 


z a,b,c,d,e,f = j- b j 

Ya,b,c,d,e,f _ Y|* b * c » d » e » f m [B]~ 1 ([A]-U) 

(16) 


2.3 FAULT SIMULATION OF A SIX PHASE POWER SYSTEM 

We are interested in calculating the fault 
generated components of voltages and currents. For 
this, the principle of superposition is used. The 
prefault values of voltages and currents can be 
found through load flow studies* Obviously, the 
fault generated voltage and current waves shall depend 
on the prefault values and the fault initiation angle* 

Consider the six phase network shown in Fig.2(a) 
Only the faulted node is shown in detail. The open 
six pole switch represents the unfaulted condition. 



^ xjO$ -< rnz mw>iT3 


S 

I 

X 


PREFAULT NETWORKS 
V PF . Rf 

^Or0r wH ^\ 




-jCbOr^S 






ab,c y de,f 


'FF 


Rg- 


(c) 


POST FAULT NETWORK 


v P a P /b ' c,d,e,f 

R f 

^Jir^ww-i 


" + Rf 


+ Rf 




f 




R g1 


(d) 


NETWORK TO EVALUATE 
FAULT GENERATED COMPONENTS 


FIG- 2 








The network remains unfaulted by connecting a six 
phase voltage source [V p p] as shown. [Vpp] represents 
the prefault voltage vector. Now the fault is 
simulated by connecting a voltage source [Vpp] with 
the polarity as shown in the figure. The value of 
[Vpp] depends on the type of fault, e.g. for a six 
phase fault [Vpp] = [Vpp]* For other types of faults 
too [Vpp] can be calculated from the fault conditions. 

To calculate only the fault generated components 
we short the voltage source [Vpp] as also all the 
independent voltage sources in the system. This is 
shown in Fig. 2(d). At the faulted bus, there is a 
series branch of the voltage source [Vpp] and the 
series fault resistance [Rp]. It is easily seen, that, 
the element of [Rp] are, 

[Rp^ij = Rf + R g i = j 

« R g i ¥ 5 

For nodal analysis, it is more convenient to convert 
the voltage source [Vpp] and impedance [Rp] to an 
equivalent current norton source. This is shown in 
Fig. 3. 

[J F ] - [Gp] [V FF ] 

[<%] « [Rp] -1 (17) 

Now the augmented bus admittance matrix for the 
power system can be derived. With [Jp] as the only 




Fig. 3 EQUIVALENT CURRENT SOURCE REPRESEN- 
TATION AT THE FAULT BUS 



f R.R - Relay point 

I 

jflG. 4 A SAMPLE 3-BUS POWER SYSTEM 

I 

i 

1 
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independent current source in the system # the fault 
generated components of voltages and currents can be 
found. It should be noted, that, all these values 
are in the transform domain and have yet to be 
converted to time domain values. 

2.4 INVERSION TO TIME DOMAIN 

Use of the Laplace domain in the past has been 
, restricted because of problems in inverting it. 

Techniques used in the past [18] have had the drawback 
of excessive storage requirements and required time 
consuming algorithms. 

A new algorithm for inverting the Laplace 
transform has been proposed based on quadrature [19]. 

We know that the Laplace transform of a time 
function u{t) is given by, 

Q0> 

F(s) = / u(t) e st dt (18) 

o 

Further, to an approximation, closed form integration 
can be replaced by a summation, l.e. , an integral 
of the type 

1 N 

/ f (r) dr = 2 W. f(r.) (19) 

-1 i=l 1 1 

It can be shown that any function f (r) in the range 
(-1, 1) can be expressed in terms of a set of Legendre 
Polynomials in this region, i.e. 
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N 

f(r) * 2 b, P k (r) r = [-l,l] (20) 

k=0 K K 

This follows from the following two properties of 
Legendre polynomials* in the region (-1,1)* 

(i) Orthogonality, and 
(ii& Completeness. 

However, by choosing and r^ suitably, we can 
make eqn, (19) to be exact for any polynomial of 
degree less than or equal to 2N-1. 

We shall show, that, 

the r^, i = 1,2. ., .N are the zeros of P N (r) 

where P^(r) is Legendre polynomial of degree N. (21) 

1 P M (r) dr 

(«>.=/ i = 1,2, .... .N (22) 

-1 (r-r^P^r^ 

To prove the first statement, consider the polynomials 
r P^j(^-) ^ = 0,1,2, .... ,N— 1. Then, 

1 , N ' 

0 « / r K P N (r) dr « 2 ®i r i P N^ r i^ k “ 0,1,2,... N-l. 

—1 i— 1 

Consider these as H linear equations in w^P^r^). Then 
the determinant, 

Jr^j i * 1,2. ...N k * 0,1,., ..N-l 

the Vandermonde matrix is non zero as all r^*s are 
distinct. Hence, the only solution is the trivial 
one, i» e. 
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® i p N tx l ) “ 0 1 “ X * N - 

Since, =9^ 0, we have, 

P N^ r i^ = 0 1 = 1 » # **** N * 

Hence, r^*s are the zeros of P^(r). 

The second equation can be verified by substituting 
it in eqn. (19). 

Finally, we make the change of variables, i.e. 

= x r = 2x - 1 

then, the quadrature formula becomes, 

1 N 

/ g(x)dx = I ok g(x ± ) g(x) = f£2x-l) 
o i=l 

x^ f s are zeros of the shifted Legendre polynomial 
P^(2x - 1 ) ^ We are interested in inverting eqn. (18). 
For this, we first map the [0, region to [0, l] 

—t 

region by using the change of variable x = e . 

Then, we have, 

F(s) *= / X s-1 g(x) dx g(x) = u(-log x) 
o 

M . 

ox 2 c^x* g(x i ) *= F(k+l) k = 0,1,2. ...N+l (23) 

v .. .... • . . '' 

Now w^’s and are known before hand for a given N. 

Hence, we can calculate these and denote them by 
matrix M. , 
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where, 

M = [w.x* ] 

1 1 ik 

We want to find g(x. ) for i = 1,2. . . .N. F(k-t-l) are 
all known for k = 0,1,2....N-1 as F(s) is known. 

Hence, 

Cg(x ± )] = M** 1 [F(k+1)] (24) 

can be calculated and stored for given N* Then, 
all that has to be done to find the inverse Laplace 
transform of F(s) is to multiply this inverse matrix 
by [F(k+1)] column vector. The inverse matrix is 
given in the Appendix along with the time instants 
tj^ * *■ log j i = To f xnd the time 

function at other time instants, use the scaling 
factor, i.e. , 

(f (at)) = Eil/ai 

2,5 NUMERICAL EXAMPLE 

A sample 3-bus power system as shown in Fig. 4 
is considered to study the fault generated travelling 
waves on a six phase transmission line* The fault is 
assumed to have occurred on the middle of one of the 
transmission lines# Identical sources are considered 
at either end of the line. 
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The one line diagram of the equivalent network 
of the sample power system is shown in Fig. 5. [A^] , 

[Bj]» [C 2 ] and [D^] are the line parameters for length 
200 km and [C 2 ] and CD 2 ] are the corres- 

ponding matrices for length 100 km. By inspection, 
the bus admittance matrix is given by 

BiVejVry -B7 1 -Bj 1 

b 7\ + b;V^s] 

-Bj 1 2B2 1 A 2 +[Y f ] 

We are interested in finding the voltage and current 
values at the relay locations. These are to be found 
in terms of [V FF ]. Solving the nodal equations, we get, 

[V 1 ] - [2B2 1 -(2B2 1 A 2 +[Y F ])B 2 (B2 1 A 1 +Bj 1 A 2 +[Y 3 ]-B" 1 ) -1 [Y F ][V FF ] 
and 

58 ®i (A r U)[V t ]. 


Y a ,b,c , d , e , f 
BUS 


-b: 


-B 


-1 


» 


Since, the travelling waves are generated only during 
the most severe faults, only the six phase to ground is 
considered. For this fault, [Vpp] = [V p p] i*e. 
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s sin 0 Q + w c cos 0 O 
s sin(0 o -6O°) + cos(0 o -6O°) 

s sin(0 o -l2O°) + u) Q cos(0 o -l2O°) 

-s sin(0 Q ) - cos 0 O 
s sin(0 o +l2O°) + {i) Q cos(0 Q +120°) 
s sin(0 Q + 60°) + (t> o cos((j) o + 60°) 

co q = natural system frequency* 

v m = maximum value of prefault voltage = 1.0 pu. 

0 O = fault initiation angle* 

Using the above [Vpp] vector, the voltage and current 
values are found at relay location 1* Then, using 
the inversion matrix, corresponding time domain values 
are found. 

The waveforms for various values of 0 Q and 
fault impedance are shown in the Figs ,6, 7 and 8* 

2.6 RESULTS AND DISCUSSIONS 

The chapter presents a frequency domain model 
suitable for study of transient phenomenon on six phase 
transmission line. The equivalent Pi-model is suitable 
for the nodal analysis of six phase power systems using 
the bus admittance matrix* The approach taken is 


[V 


FF- 


V m 

ZTS 
S +oj_ 
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age waveforms at relay 



Fig. 7 Fault generated voltage waveform at relay 
location (0 = 90°) 





xjot- 4 xioot ^ 


Fig. 8 Fault generated voltage waveform, at relay 
location (0 = 0°) 
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completely general and can be used for study of faults 
on large systems too. 

The chapter presents a novel technique for 

inversion of the Laplace transform. The technique 

based on quadrature principle shows hew the inversion 

process can be approximated as a matrix multiplication. 

The number of operations required to find the time 

2 

function at n time instants varies as n . It is clear, 
that, the technique is superior in many respects to the 
conventional Fourier transform techniques. Further, 
no complex numbers are involved thus saving computa- 
tional time further* Yet another advantage of the 
scheme is, that, it does not require much storage 

space. For example, the inverse matrix requires only 
2 

the n matrix elements to be stored along with the 
n zeros of the shifted Legendre polynomials. It is 
hoped, that, this technique would be found useful in 
transient analysis of power systems* 



CHAPTER 3 


PROPOSED UHS RELAYING SCHEME FOR SIX PHASE TRANSMISSION 

LINE BASED UPON TRAVELLING WAVE PHENOMENA 

3.1 INTRODUCTION 

Designing six phase systems would require 
alongside, adequate protection schemes which are 
fast in operation and reliable. Over the years, 
various protection schemes have been proposed based 
on steady state analysis as well as transient 
behaviour. In the recent years, travelling wave 
relaying schemes have been proposed. These relaying 
schemes come in the category of ultra high speed 
schemes, i.e* (it operates in less than quarter cycle). 
This ultra high speed clearing of fault would have a 
two fold advantage i) it would improve the 
transient stability of the system as the kinetic 
energy introduced into the system is proportion to 
the square of the fault clearing time [20], 

Secondly, this also improves system reliability. 

In the past, the feasibility of ultra high speed 
relaying scheme was not very bright because of the 
absdnce of high speed circuit breakers. However, 
with recent advances in circuit breaker technology 
[1,2] attention is now focussed on ultra high speed 
protective relays. 
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These ultra high speed relays are based on the 
travelling wave phenomena which are generated on the 
occurrence of a fault and travel towards relays at 
either end from fault location* 

Takagi et al. [10] and Akimoto et al. [11], 
defined a relaying signal, using sending and 
receiving end voltages and currents and also the 
surge impedance of the line* The main demerit is 
that the scheme requires exchange of information 
between ends of the protected line. Fault location 
* scheme proposed by above authors are to complicated 
for easy realisation* 

An ultra high speed relaying scheme developed 
by ASEA, Sweden has been tested on Bonville Power 
Administration (U.S,A. ) (BPA) 500 kV system* It needs 
the exchange of only qualitative information between 
the ends of the protected line. Basically, the scheme 
is a directional comparison scheme. Fault location 
schemes have been proposed by Vittins [6], Fault 
direction determination algorithms have also been 
proposed by Johns [9] and Agarwal* 

In this chapter, an amplitude comparison 
relaying scheme based on travelling wave phenomenon 
has been proposed. The scheme utilises the reflection 
properties of travelling waves to derive appropriate 
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relaying signals [21] . The scheme proposed is simple 
and reliable, A major advantage is, that, it does not 
require any transfer of information between the ends 
of the protected line. The results of digital 
simulation of a faulted sample power system are 
presented to confirm the feasibility of the proposed 
scheme, 

3.2 PRINCIPLE OF OPERATION 

The principle of operation can be best studied 
by considering an equivalent single phase system. 
Consider the single phase transmission line connecting 
two large power systems as shown in Fig. 9* After 
the occurrence of the fault, the voltage and current 
at any point in the system can be regarded as the sum 
of prefault and fault generated values, i.e. 


v 

i 


v + v f 

v + % 


The fault generated components, travel from the fault 
location towards either end of the protected line. 
Assuming a loss less line, the fault generated 
components can be represented as 


v f * f x (x-ut) + f 2 (x+ut) 
i f * 2T * f 2 (x-ut)) 
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S-E. R.E. 



FIG- 9 A TRANSMISSION LINE INTERCONNECTING 
TWO POWER SYSTEMS 


S-E. F R.E. 



FIGio LATTICE DIAGRAM FOR AN INTERNAL 
FAULT 
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where , 

u = velocity of propagation on the line 
Z Q = surge impedance of the line. 

These travelling waves on reaching the relay 
locations undergo reflection and refraction. If p s 
and be the reflection coefficients at either end 
of the line, the subsequent waves on the line can be 
represented by a Bewley diagram shown in Fig. 10. v ff 
is the fault generated component of the voltage at 
the fault point* For a bolted fault it is the prefault 
voltage at the point of fault. 

The proposed relaying scheme uses the following 
pair of signals at each relay location. For the 
sending end relay, the input signals in p*u. would be 


S 1 

" v fs “ R set 

^■f s 

s 2 

“ v fs + R set 

Hs 


v^ s and i|. g are the fault generated components of 
voltage and current at sending end, and, 

R. .. = surge impedance of the line. 

m W m 

For the faulted line, the values of v fg , i fg , JSjJ 
and jS 2 | at various instants have been shown in Table 1* tj 
represents the time of travel of the wave fr am the fault 
location to the relay location. Since magnitude of p s 
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is less than 1.0, it is seen that Js^ j ) |S 2 | fro® 
t = T i till a further reflected wave appears at the 
relay location. 

Now, consider a fault which occurs outside 
the protected zone. When the travelling wave schemes 
the relay location, it gets reflected and refracted 
at the sending end bus* Due to the refracted wave 
we have, 

v fs “ f t v ff ^"f s = j?t v ff/ Z o 

where 

p ^ « refraction coefficient at the sending end 

bus 

** surge impedance of the faulted line. 

Then, we have, 

s i - 

Z o 

S 2 “ ^ + V ff (27) 

o 

Hence, for some duration, in this case, jS^J < |S 2 i » 
Thus we have the following result, 

JSjJ > |S 2 j for internal fault 

|S^ J < )S 2 j for external fault. (2.8) 

This property is used for designing an appropriate 
relaying scheme. 
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3.3 EXTENSION TO A SIX PHASE SYSTEM 


From equations (5,6) we know that the transmission 
line equation of the six phase system can be written 
in decoupled form using the symmetrical transformation 
matrix [T], i.e. 


- oV (m) 
dX 2 m 

m 0 y I f 2 t 3 f 4 t 5 


d 2 I (m) = a !(«>) 

o m 


(29) 


dX' 


This suggests, that, the modes propagate independently 
of each other and that a six phase transmission line 
behaves as a single phase line for each of these 
modes separately. It follows, therefore, that the 
results obtained in the previous section, apply 
equally well to each modal component separately. 

Hence, six pairs of signals, one for each mode, are 
employed. For the sending end relay, these are given 


by. 



,(ra) 

r fs 


R (m) Am) 

H set is 

hi 3=5 0 i l i 2 i 3 i 4 f 5 


4 m) ~ + < 3 °> 

where, v^®^, qre the modal component values of 

fault generated voltages and currents. These are 



RELAY SIGNAL DIFFERENCE ( p.u.) 



FIG. ii Relay signal difference iSil- IS2I for internal 
fault on sample network. 
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found using the symmetrical component transformation 

matrix [T], 

(m) 

R set is the Sur 9 e impedance for the m th mode. 

The amplitude comparison scheme is applied to each 
signal pair. The tripping or blocking signal is given 
by the first positive outcome of the comparison, i.e. 

Js[ m) | > \S^\ m = 0,1,2, 3,4,5 ! INTERNAL FAULT 

js£ m) | < |s| m > j : EXTERNAL FAULT 

(31) 

3.4 TESTING ON A SAMPLE POWER SYSTEM 

The above scheme is tested by simulating fault 
conditions on the sample powers system of Chapter 2, 

The signals and S 2 are shown in Figs, 11,12 for 
internal and external faults. The values of signal 
difference confirm the above principles, 

3.5 IMPLEMENTATION OF THE RELAYING SCHEME 

With the emergence of fast and reliable 
microprocessors and dedicated software, digital 
protection of transmission lines has become feasible. 
Digital protection of fears the two fold advantage of 
speed and flexibility. 

Digital travelling wave algorithms use the 
post fault transient current and voltage waveforms 
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which are obtained after filtering out the prefault 
sinusoidal voltage and current using a simple digital 
filter. The relaying decisions are thus based on the 
fault generated components of voltage and current. 

Very few algorithms based on travelling wave phenomenon 
have been proposed. This section suggests the algorithm 
for implementing the relaying scheme proposed in the 
earlier sections, 

3*5.1 Fault Detection Algorithm 

This is similar to the algorithm proposed by 
Mann and Morrison [12 ], One counter is provided for 
each phase voltage* The difference between currently 
sampled instantaneous voltage of each phase and the 
corresponding value of the previous cycle is stored. 

If this difference exceeds a pre set value (0.05 pu) 
in any case the counter of that phase is incremented 
by one. If not,, the counter is decremented by one. 

Once the counter reaches a pre set value in any of 
the six phases, fault is assumed to have occurred and 
the relaying algorithm is executed. 

3.5.2 Relaying Algorithm 

Once the fault detection algorithm yields a 
picture outcome, the relaying algorithm Is executed. 

This is as follows 
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Step 1; 


Step 2 s 


Step 3: 
Step 4: 


Step 5s 


Step 6s 


Calculate by cycle to cycle comparison, the 
fault generated components of all line 
currents of five sample sets backwards from 
the set at which fault detection algorithm 
gave a logical yes output. 

Calculate the modal components of phase 
voltages and line currents. At the first 
instance, the first sample set is considered. 
Execute steps 4 to 5 for the modes 0,1,2, 3, 4,5. 
Check if lsj m ^j > js| m ^J. If so, start 
the transmitter. If trip signal arrives 
from other end, give the trip signal to 
circuit breaker. 

If JS^ m > | < |s| m) j , block the transmitter 
signal. 

Take the next sample set and go to step 2. 


3.6 DISCUSSIONS AND CONCLUSION 

The chapter presents a new relaying scheme for 
protection of multiphase transmission lines. The 
schemes use the presence of travelling wave component 
in the post fault waveforms to derive the appropriate 
relay input signals. The proposed scheme is simple 
and does not require complex computations. 
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The viability of the scheme has been demonstrated 
on a sample 3— bus power system network. The chapter 
also presents an algorithm for digital implementation 
of the scheme. 

The relaying schemes based on travelling waves 
are immune to power swings since the latter are 
comparatively slow transients. The transients due to 
switching appears as an external disturbance and hence 
will not operate the relay. Since the fault detection 
takes place before the CT*s saturate, the associated 
problems are eliminated. The lightning strokes have 
a high frequency component associated with them in 
the MHz range# These could be eliminated using a 
simple lowpass filter. 

The fault detection algorithm proposed in this 
chapter can be used also for various analog travelling 
wave relaying algorithms for digital implementation. 
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TABLE 1 


Relay Input Signals for an Internal Fault 
Time v f 8 if S jS 2 j 

t <t x 0 0 0 0 

* = "C v ff * v ff/^o ^ l v ff i ^ 

4 >^1 {1+ )s\ff ~ (1 *fs )v ff / ' Z o ) 2 l v ffl 2 lfs v ffl 

till another 
wave pair 
arrives 


CHAPTER 4 


CONCLUSION 

The ever increasing demand for power have 
prompted power engineers to search for newer and 
more efficient methods of power generation and 
transmission. In this context, high phase order 
transmission (HPOT) system seems to be a promising 
way of transmitting greater amounts of power to the 
load centres at reasonable cost and at better 
voltage profile. 

Any proposalfor design of a multiphase power 
system must be accompanied by adequate measures 
against faults. The ultra high speed relaying 
schemes based on travelling wave phenomenon offer 
the dual advantage of speed and reliability. Hence, 

IMS protection is particularly suited to protection of 
EHV/UHV transmission lines. The use of such schemes is 
particularly facilitated by the advent of high speed, 
low cost microcomputers and microprocessors. 

Accordingly, the thesis aims primarily at 
designing protective relaying schemes for six phase 
transmission lines based on travelling wave phenomenon* 
Fox this a frequency domain Pi-model of a six phase 
transmission line has been developed. The thesis has 
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also explored a novel technique for Laplace transform 
inversion which, it is hoped, shall facilitate 
transient studies on power systems. 

Designing a relaying scheme for a multiphase 
transmission line requires that a model for the trans- 
mission line which would lend itself to easy transient 
analysis be developed first. To this end, a frequency 
(S-domain) model for the transmission line has been 
developed. The model has been developed keeping in 
mind, that it should be helpful in studying generation 
and propagation of travelling waves on the transmission 
line. Laplace domain techniques have been used for 
this purpose. 

In the past, because of difficulties in inverting 
a Laplace domain function to the time domain, use of 
Laplace transform was somewhat restricted. The thesis 
has suggested a new method of inverting Laplace 
transforms. The technique does not require any complex 
operations like contour integration etc. and simply 
reduces the Laplace inversion to a matrix multiplication. 
The matrix in question, is a constant, and can be stored 
for repeated use. The thesis also presents a generalised 
nodal analysis for fault analysis of multiphase power 
systems. The S— domain bus admittance matrix has been 
used for the purpose. Similarly, the Pi-model of the 
transmission line is also suited to nodal analysis. 
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Th@ viability of model developed has been 
demonstrated on a sample 3-bus system. The voltage 
waveforms for different values of fault resistance 
and fault Initiation angle 0 have been shown 

high speed clearing of faults improves 
the transient stability of the systems. The fault 
clearance time depends on both the relay speed as well 
as delay associated with the circuit breaker. With 
development of ultra high speed circuit breaker, ultra 
high speed relaying schemes have gained importance. 

But very few relaying schemes based on travelling wave 
phenomenon has been proposed so far. This paper 
present® an amplitude comparison relaying scheme based 
on travelling wave phenomenon along with fault detection 
algorithm. 

The amplitude comparison scheme decomposes the 
fault generated voltages and currents to their modal 
components and makes trip/block decisions based on 
these modal components. Tripping is initiated if the 
fault is detected as a forward fault at both ends. 

The underlying principles have been verified by fault 
simulation on a sample power system. The effect of 
fault resistance and fault initiation angle has been 
included in the study. With the use of modal components 
all types of faults can be protected against, using the 
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same algorithm* Hence the relaying scheme can be 
considered to be a new type of polyphase relay. 

Digital protection of transmission line has 
largely been restricted to various types of distance 
relay algorithms which involve the computation of 
power frequency Impedance from the complex post fault 
waveforms of voltages and currents. Different types 
of analog and digital filters have been constructed for 
extract ingthe fundamental frequency components of 
voltages and currents. However, the accuracy of 
determining this component depends on the extent of 
filtering. Further, the nature of filtering depends 
upon the rigour of modelling the line. Either more 
elaborate modelling or complex filtering can result 
in inordinate time delays. But travelling wave relays 
require that only the prefault component of voltages 
and currents be filtered out. Further, only a few 
samples are required for this purpose. Hence, the 
speed of operation of UHS relays is much faster. 

The prefault voltages have been filtered out 
using the cycle to cycle comparison method proposed 
by Mann and Morrison. The viability of these algorithms 
has been tested on Dec 1090. 

Future Scope of Work i The modelling of transmission 
line, fault simulation and testing have been done using 
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software techniques. Such techniques, while demonstrating 
the theoretical feasibility of schemes, do not bring 
forward the practical difficulties encountered during 
hardware, on-line implementation. 

Moreover, the thesis has considered only simple 
multiphase systems. Future work could include the 
modelling of transformers, voltage and current transducers , 
generators in the overall system to obtain a more realistic 
picture. 

Lastly, relaying principles developed shall need 
reassessment when applied to double circuit lines, series 
and shunt compensated lines and multi-terminal lines. 
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APPENDIX - A 


THE ELEMENTS OF THE INVERSE MATRIX WITH DIVISION 

BY WEIGHTS 


In these tables, we present the inverse of the 

matrix (w^x^* 1 ), i,j = i, 2,.».,,N f for N = 6. 

N = 6 

ROW 1 1.827E 1 

— 2.263 E 2 
9.730E 2 

— 1.886 E 3 
1.691E 3 

-5.702E 2 

ROW 2 -5.213E 0 

1.881E 2 

-1.078E 3 

2.390E 3 

-2.309E 3 

8.160E 2 

ROW 3 2.636E 0 

-1.038E 2 

8.347E 2 

-2.237E 3 

2.429E 3 

— 9,275 E 2 

ROW 4 -1.620E O 

6.546E 1 

-5.750E 2 
1.794E 3 

— 2.208E 3 

9.275E 2 

ROW 5 1.063E O 

— 4.337 E 1 
3.943E 2 

-1.311E 3 

1.770E 3 

-8.160E 2 

-6.387E -1 
2.616E 1 

-2.411E 2 

8.234E 2 

-1.159E 3 

5.702E 2 


ROW 6 
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APPENDIX - B 

NEGATIVES OF LOGARITHMS OF ROOTS OF SHIFTED 
LEGENDRE POLYNOMIALS 

The negative logarithms tabulated are 

t ± = -In x i 

where 

P N (Xi) = 0 for i = 1,2,.. ,,N} N = 6 

N = 6 
3*388323 
1.775520 
0*965769 
0.479150 
0.185601 


0.034348 
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